Identifying Dsj(2700) through its decay modes 
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We study how to assign the recently observed Dsj (2700) meson to an appropriate level of the cs 
spectrum by the analysis of its decay modes in final states comprising a light pseudoscalar meson. 
We use an effective lagrangian approach with heavy quark and chiral symmetries, obtaining that the 
measurement of the D* K decay width would allow to distinguish between two possible assignments. 

PACS numbers: 13.25.Ft,12.39.Fe,12.39.Hg 



I. INTRODUCTION 

Recently, Belle Collaboration collected an important 
new piece of information on charm spectroscopy from 
the first Dalitz plot analysis of the decay process £?+ 
D°D°K+ The Dahtz plot shows an accumulation of 
events at M^{D°D°) = 16-18 GeV^ and AP{D°K+) = 
7 — 8 GeV^. This could be interpreted as the overlap 
of two contributions, a (horizontal) band in the D^D^ 
channel due to a state possibly identified with -0(4160), 
and a (vertical) band in the D^K^ channel due to a cs 
state which does not correspond to any of the previously 
observed mesons with this quark content. A horizontal 
band at AP{D°D°) ~ 14.2 GeV^ is also present, due to 
the production of V'(3770). 

Information from the Dalitz plot was enriched by the 
study of the individual invariant mass distributions. The 
AP{D'^K^) distribution shows the presence of a reso- 
nance, Dsj{2700), with parameters: 

M = 2708 ± 9tlo MeV 
r = 108±23t^'i' MeV. (1.1) 

Moreover, from the distribution in the helicity angle 6, 
the angle between the momentum and the opposite 
of the kaon momentum in the D'^D'^ rest frame, it was 
possible to assign the spin-parity = 1^ to Dsj (2700). 
A measurement of the branching fractions 

B{B+ ^ D°D''K+) = (22.2 ± 2.2Uattl%ys) x 10"^ 
8(3+ (2700)5°) X ^(£1^7(2700) D°K+) = 



{n.3±2.2\statt2.t\sys) X 10"^ 



(1.2) 



was also carried out. 

Belle's observation is just the last one of a scries of dis- 
coveries of open and hidden charm hadrons which have 
greatly enriched the charm spectroscopy in the past few 
years 0]. In case of open charm, which _Dsj(2700) be- 
longs to, before the B-factories era the known cs spec- 
trum consisted of four states only: the pseudoscalar 



£'s(1968) and vector £'*(2112) mesons, corresponding 
to the s-wave states of the quark model, and the axial- 
vector Dsi{2536) and tensor Ds2{2573) mesons, p-wave 
states. In 2003, two narrow resonances: Z?sj (2317) and 
Dlj{2A60) were discovered by BaBar [3| and CLEO [i| 
Collaborations, respectively, and later on confirmed by 
other experiments [5|. They were assigned spin-parities 
= O"*" and = 1"''. Their identification as proper cs 
states has been controversial since then [6|; however, they 
have the quantum numbers of the two states needed to 
complete the p-wave multiplet, and their radiative decays 
occur accordingly, so that their interpretation as ordinary 
cs configurations is natural 0- 

Last year, BaBar Collaboration announced the ob- 
servation of another cs meson, Dgj (2860) decaying to 
D°K+ and D+Ks, with mass M{Dsj{2860)) = 2856.6± 
1.5 ± 5.0 MeV and width r(L),j(2860)) = 47 ± 7 ± 10 
MeV Q . The quantum numbers of this meson have still 
to be assigned, and two proposals have been put for- 
ward. One is that this state is the first radial excitation 
of Dsji2317), hence with = 0+ the other one 

is that the meson has = 3~, an interpretation 
posed in JTj] and supported by a lattice QCD study 
In the same analysis of the DK mass distribution, BaBar 
noticed a broad structure with M = 2688 ± 4 ± 3 MeV 
and r = 112±7±36 MeV [8] , likely the same resonance 
£13 7(2700) found by Belle. 

In this scenario, Dsj (2700) needs to be properly placed 
in the spectrum of mesons with charm and strangeness. 
In our study we investigate its decay modes with the aim 
of finding signatures allowing to assign the meson to a 
particular cs level. In this respect we follow a different 
strategy from the one adopted in the framework of quark 
models flo| , where the assignement is made by comparing 
the observed mass to the predicted value in a suitably 
chosen interquark potential, a procedure which can find 
difficulties in the cases where neglected effects, such as 
the threshold effects, are relevant. 

In Section |TT] we present a theoretical framework based 
upon an effective Lagrangian describing strong decays 
of heavy mesons to final states comprising a light pseu- 
doscalar meson, and displaying heavy quark and chiral 
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symmetries. We analyze the various decay modes and 
compare the predictions that follow from different assign- 
ments. The implications for the still unobserved states, 
including cq mesons, are discussed in Sections IllUfVl then 
our conclusions are presented. 



II. HEAVY MESONS DECAYS TO LIGHT 
PSEUDOSCALAR MESONS 

The study of properties and decays of hadrons con- 
taining a single heavy quark Q — c, b is suitably carried 
out in the mg —^ oo limit, which is formulated in the 
Heavy Quark Effective Theory [l^. In such a limit, the 
heavy quark acts as a static colour source for the rest 
of the hadron; its spin sq is decoupled from the total 
angular momentum of the light degrees of freedom s^, 
and they are separately conserved. Hadrons can be clas- 
sified according to the values of S£ and of the total spin 
J = Sq + S£. In particular, heavy mesons can be orga- 
nized in doublets, each one corresponding to a particular 
value of sg and parity; the members of each doublet differ 
for the orientation of sq with respect to Si and, in the 
heavy quark limit, are degenerate. Mass degeneracy is 
broken at order l/mg. 

For Qq states one can write sg — Sq + £, where Sq 
is the light antiquark spin and i is the orbital angular 
momentum of the light degrees of freedom relative to the 
heavy quark. The lowest lying Qq mesons correspond 
to £ = (s-wave states of the quark model) with sf = 

^ . This doublet comprises two states with spin-parity 
= (0^,1^). For £ = 1 (p-wave states of the quark 
model), it could be either sf = \ or = | , the 



i or sf = I 
two corresponding doublets having ~ (0+, 1+) and 
= (l+,2+). The mesons with £ = 2 (d-wave states) 
are collected either in the sf ~ ^ doublet, consisting 
of states with = (1^,2^), or in the sf = ^ doublet 
with = (2^,3^) states. And so on. 

The two states i:)s(1968) and D*{2n2) can be identi- 
fied with the members of the lowest lying sf — ^ dou- 
blet. The resonances iDsi(2536) and Z)s2(2573), together 
with Dsj{2317) and 15*^(2460), fill the four p-wave lev- 
els: in particular, Z?s2(2573) corresponds to sf = 

JP = 2+ state, while Dsj(23l7) to sf = i^, = 0+. 
The two axial- vector mesons Dsi{2536) and Z?*j(2460) 
could be assigned to a linear combination of si' — |^ 

and sf — states, which is allowed at ©(l/mg); how- 
ever, in case of non-strange charm mesons such a mixing 
has been measured and found to be small [3, [iBl, so 
that we can identify Dsi{2536) and £);j(2460) with the 



JP = 1+ sf = i+ and sf 



1 + 



2 ^ii^ — 2 states, respectively. 
In the interpretation in Dsj{2860) corresponds to 
the jP = 3^ component of the sf — ^ doublet, so that 
the still unassigned levels are the two states belonging 
to the = I doublet, with = (1~,2~), and the 



partner of Dsj (2860) in the sf = ^ doublet, which has 
spin two. 

The classification can be continued for higher values 
of sf to describe high spin mesons; moreover, it can be 
replicated for radial excitations, and it is expected to hold 
as far as ©(l/mg) effects are small. 

Since the spin-parity of Z?sj (2700) has been deter- 
mined: = 1~, the state fits either in the doublet 
with sf = ^ or in the one with sf ~ ^ . However, 

a 1^ state belonging to the sf = ^ doublet is already 
known, D*{2112), so that in this case Dsj(2700) would 
be a radial excitation. Combinations of the two cases 
are not allowed in the heavy quark limit; a possible role 
of 1/mg effects will be discussed below. Therefore, two 
possibilities must be considered: 

• Dsj(2700) belongs to the doublet with sf ^ ^ 
and is the first radial excitation; we denote this 
state as D!'; 



• (2700) is the low lying state with sf 
denoted as D*^. 



These two possibilities can be conveniently analyzed 
adopting a formalism which represents the various dou- 
blets by effective fields being 4x4 matrices [l^ . The two 



doublets of interest, with si 



denoted by Ha, and 



with sf — ^ denoted by Xa (a is a light flavour index) 
are given by: 



Ha 



(2.1) 



with V the meson four-velocity. Notice that Ha describes 
the fundamental sf — ^ doublet; the doublet corre- 
sponding to the first radial excitations is described by 
an identical structure H'a- The various operators in (|2.1|) 
annihilate mesons of four-velocity v which is conserved 
in strong interaction processes: the heavy field operators 
contain a factor y/mp and have dimension 3/2. 

The interaction of these heavy mesons with the octet 
of light pseudoscalar mesons, introduced using the fields 
^ — e , with the matrix A4 containing tt, K and rj fields: 



M = 



V 



TT 

K- 



K+ 



(2.2) 

and /tt = 132 MeV, can be described by an effective 
Lagrangian invariant under chiral transformations of the 
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light fields and heavy-quark spin-fiavour transformations 
of the heavy fields ■ At the leading order in the heavy 
quark mass and light meson momentum expansion the 
decays F HM [F = H' and X, and M a light pseu- 
doscalar meson) can be described by the Lagrangian in- 
teraction terms [l7j : 

Ch' = ~gTr[H',H,^,j,AL] (2.3) 
Cx = ^Tr[HaX^{iD^A + i P A^^^alz] + h.c. (2.4) 

where Af^ba = § ("^^^^^ - C^^^Oba ^'^'^ ^ 
variant derivative D^^a = -had^j, + V^ba with V^ta = 
\ (^^i9^C + ?^M^^)ba- chiral symmetry-breaking 

scale and g , k' are effective couplings. We use = 1 
GeV, while at present the value of the couplings is not 
known neither from experiment nor from theoretical con- 
siderations. 

The widths of the allowed decay modes of Dsj(2700) 
for the two possible interpretations D*J or D*^ can be 
computed using (|2.3p and (|2.4p . Since a heavy vector 
state can decay to a light pseudoscalar and a heavy 
pseudoscalar or vector meson, we consider the modes: 
D^j(2700) ^ D+Ks, D^K+, D,r] and D,j(2700) ^ 
D*+Ks, D*°K+, D*ri, and compute the ratios of decay 
widths: 



i [Us J UJ\) 

T{Dsj ^ Ds-n) 

= nDsj ^ DK) ^'-'^ 

r(n. , n*'n^ 

R3 





^ D*K) 




DK) 


r{Dsj - 




TiDsj - 


^ DK) 


r(A./ - 


- d:v) 


TiDsj - 


^ DK) 



(with D'^*^K = D^*'>+Ks + D^*'>°K+) for both the as- 
signments of sf to Dsj{270Q). These ratios are useful 
to discriminate between the two assignments; moreover, 
the method is sensible, since in the ratios (|2.5p the de- 
pendence on the (unknown) effective couplings drops out, 
and the predictions are model independent. 

Let us first identify Dsj{2700) with £)*'. Using the 
effective Lagrangian (|2.3p . the I?*' decay widths can be 
written as follows: 



T{D:' ^ DaPa) 

t{d:' ^ d:p,) 



^2 



Ci 



M 



Da 



67r/2 Mzj. 
37r/2 Md^ 



(2.6) 



where the light flavour index a = u,d,s identifies 
Di*) = £)(*)+, Di*^ respectively, and Pa repre- 

sents a light pseudoscalar meson with quark content sa 
{Pa = K^, Ks{L), v)- Cp is a coefficient depending 

1 2 

on the P meson: Cx+ = 1, Cks — 7: and = -; 

^ o 

the modulus of the three momentum q reads: \q\ = 
Ai/2(M2.„M2,.,,M|,J/2Mz3j.. 



On the other hand, if Dsj{27QQ) is identified with D*^^, 
from the Lagrangian (|2.4|) we obtain: 



r(Z?:i ^ DaPa) = 



C. 



16 



k' 



97r/2 \A^J Md: 
T{D:, ^ DlPa) 



^ [Ml + kf ] |gt 



k' 



97r/2 \K^) Md.^ 



+ kf]kf- (2.7) 



These two sets of expressions produce different values 
for the ratios (j2.5p . as one can appreciate considering 
Table |T] where the numerical results are collected (with 
the errors obtained considering the uncertainty in the 
Dsj(2700) mass). 



TABLE I: Ratios Ri for Dsj{27QQ) identified as D*' or Z)*i. 





Ri X 10^ 


R2 X 10^ 


R3 X 10^ 




91 ±4 


20 ± 1 


5 ± 2 




4.3 ±0.2 


16.3 ±0.9 


0.18 ±0.07 



The ratios Ri and R3 are very different if -Dsj(2700) 
is I?*' or D*i, so that the measurements of these ratios 
allow to properly identify the Z?sj(2700). Since R3 is 
small, the D* K decay mode is the main signal that must 
be investigated in order to distinguish between the two 
possible assignments for Z?sj(2700). 

It is worth noticing that the results in Table |T] are dif- 
ferent from the ones obtained in ref. jl^] using the ^Pq 
model, a quark model with harmonic oscillator meson 
wave functions, where it turns out that the DK mode is 
suppressed if -Dsj(2700) is identified with D*J . Moreover, 
in the D^^ identification, more reduction of the D*K sig- 
nal than obtained in ref.[l^ is expected. 

From the ratios in Table IH assuming that the width 
of Dsj{27QQ) is saturated by decay modes with a heavy 
meson and a light pseudoscalar meson in the final state, 
we can determine the coupling constants appearing in 
(|2.3p and (|2.4p , in correspondence of the two possible as- 
signments for Z3sj(2700). This assumption is reasonable, 
since decay modes with more than one light pseudoscalar 
meson, or a light vector meson in the final states are 
severely phase-space suppressed. 

Identifying £'^,7(2700) with D*J we obtain: 



g = 0.26 ±0.05 
while if Dsj(2700) is D*^^ we obtain 
k' = 0.14 ±0.03. 



(2.8) 
(2.9) 



These two values are similar to the results obtained for 
analogous coupling constants appearing in the effective 
heavy quark chiral Lagrangians [23|. Using the results 
(I2.8p - (|2.9p the branching fractions of the various decay 
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modes can be computed; the results are collected in Table 
im The errors in eqs. (|2.8p and (|2.9p . as well as those in 
Table are obtained from the uncertainties in , 
and r{Dsj). 

A comment is in order about the accuracy of the re- 
sults. The effective Lagrangians (|2.3|) - ()2.4p coincide with 
the first terms of an expansion in the light pseudoscalar 
meson momenta. Since in the decays we have considered, 
such momenta are not very small, one should in princi- 
ple add other terms, which should be weighted by new, 
unknown cupling constants. Analogusly, corrections to 
the heavy quark hmit , which is also used in (|2.3p - (|2.4p . 
could be considered, by adding 0(\/mQ) terms, which 
would contain new unknown constants as well [la]. We 
cannot assess the role and the size of such corrections 
on general grounds, however we expect that they would 
largely cancel out in the ratios of widths. On the other 
hand, a mixing between the two — 1^ , = ^ , | 
states, which is possible at 0{l/mQ), would involve a 
mixing angle that could be fixed by measuring the ratios 
in Table |T] and comparing the experimental results with 
the ratios computed in the heavy quark limit. 



III. PARTNERS WITHOUT STRANGENESS 

It is interesting to study the charmed mesons with the 
same quantum numbers as Dgj (2700), but with a dif- 
ferent light quark flavour. These states are a charged 
charmed meson and a neutral one, denoted as and 
£)j, respectively. They have not been observed yet, so 
that their masses are unknown. We fix such masses 
to 2600 ± 50 MeV by the reasonable assumption that 
Z)sj(2700) is heavier by an amount of the size of the 
strange quark mass. 

Allowed decay modes for L»|(2600) are: D] 0"^+, 
D+TT°, DsK°s(L), D+T], and Dj D*"tt+, D*+tt°, 
D*+-q, while for D^j they are: D^, D+n-, D°tt°, 
D,K-, D^rj and 0° D*^7r°, D*+n-, D*°rj; the cor- 
responding widths are obtained using eq. (|2.6ll2.7p . and 



+ (0) 



The states 



depend on the possible identification of D 
having — 2 ^'^'^ denoted as and are radial ex- 



citations, while the states having si 

*+(o) 



are denoted 



as D 

Using the effective coupling constants g and k' in 
^M,^M, we obtain: 



r(£)*'+(o)) ^ (128 ± 61) MeV 



T{Dl^^°'^) = (85 ±46) MeV 



(3.1) 



(3.2) 



so that the cq partners have widths wich are different in 
the case of the two assignments, although with a size- 
able uncertainty. Since the mesons are not very broad. 



it should be possible to observe them. The predicted 
branching fractions, collected in Table Hill confirm that 
the two assignments produce different results. In the 
identification with the state _D*', the mode D*' D*tt 
has the largest branching fraction, while in the second 
hypothesis, i.e. Dj = D'l, the mode with the largest 
branching ratio is Dl — )■ Dn. 



IV. SPIN PARTNERS 

Since in the heavy quark limit the heavy mesons are 
collected in doublets with a definite value of sf , the state 
Ds J (2700) has a partner from which it differs only for the 
value of the total spin. 

The partner of D*' {sf = - ) has = 0" 



the first radial excitation of D, 

3". 



it is 
5. On the 

has = 



denoted 

other hand, the partner of D*i {sf 

2^: we refer to this state as to D*2. In both cases, the 
decay modes D'^, D*^ D*°K+, D*+K°^^y D^r], are 

permitted. Using the effective Lagrangians 
find: 



we 



m's- 
m*s2 ' 



-9" 



Mr 



2^/2 Md'^ 



A. 



D' 



Mi 



(4.1) 



In the heavy quark limit, these partners are degener- 
ate, hence, in the numerical analysis we assign them the 
same mass as Dsj (2700). Using the obtained values for 
g and fc', we get: 



no's) 

r(^:2) 



(70 ± 30) MeV 
(12 ± 5) MeV 



(4.2) 



and the branching fractions in Table [TVl In the two as- 
signments the spin partners differ for their decay width. 



V. J (2700) DECAY CONSTANT 

We conclude our study observing that, since -Dsj(2700) 
has been discovered through the production in 5+ 

Ds j{2700), it is possible to estimate its decay constant 
Jd^j defined as: 

< 0|s7^(l - j5)c\Dsjip, e) >- fD^,, Md,,, . (5.1) 

The effective Hamiltonian governing b — )■ ccs transi- 
tions: 

H^ff = ^K^K; [CMOi + C2{^i)02] , (5.2) 
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TABLE 11: Ds J {2700) branching fractions in correspondence to the two assignments. 





B{Dsj D°K+) 


B{Ds., ^ D+Ks) 


B{Ds.j ~* DsTj) 


B{Dsj D*°K+) 


B{D,j ^ D'+Ks) 


B{Dsj -> D^T^) 




(24 ± 14)% 


(12 ±7.0)% 


(7 ±4)% 


(22 ± 13)% 


(10 ± 6)% 


(1.7 ±1.2)% 


Dti 


(44 ± 25)% 


(21 ± 12)% 


(11 ±6)% 


(1.9 ±1.1)% 


(0.9 ± 0.5)% 


(0.12 ±0.09)% 



TABLE IIL Branching ratios of the cq partners of Dsj (2700), and Dj, for the two possible assignments; the mass of the 
states is fixed to 2600 ± 50 MeV. 





B{Dj ^ L»%+) 


B{D+ D+tt") 


B{Dj ^ DsK^s) 


B(D+ ^ 


B{D+ D*\+) 


B{D+ D'+n°) 


B{D+ ^ D*+rj) 


D*'+ 


(27.0 ± 2.1)% 


(13.3 ± 1.0)% 


(2.3 ±0.8)% 


(5.3 ±1.0)% 


(32.4 ± 0.8)% 


(16.1 ±0.4)% 


(1.2 ±1.8)% 




(51.1 ±3.5)% 


(25.0 ± 1.7)% 


(4.0 ± 1.4)% 


(11.8 ±2.1)% 


(2.7 ±0.1)% 


(1.3 ±0.1)% 


(0.1 ±0.2)% 




B{D°, ^ D+-K-) 


B{D°j D"-k") 


B(I>"j ^ DsK-) 


B{D'} ^ D\) 


B{D"j D'+n-) 


B{D'} ^ D*V) 


S(Z)5 ^ D*°v) 




(26.5 ±2.1)% 


(13.5 ± 1.1)% 


(4.9 ± 1.6)% 


(5.5 ±1.0)% 


(32.0 ± 0.7)% 


(16.3 ±0.4)% 


(1.3 ±1.9)% 


Df 


(49.7 ± 3.3)% 


(25.6 ± 1.8)% 


(8.3 ±2.8)% 


(12.3 ± 2.0)% 


(2.6 ±0.1)% 


(1.3 ±0.1)% 


(0.1 ±0.2)% 



where C'l and C'2 are Wilson coefficients and penguin 
operators have been neglected, involves the operators: 



p2 = 1.179 ± 0.048 ± 0.028 from a recent BaBar measure- 



01 = [C74l-75)&][s7''(l-75)c] 

02 = [c7^(l-75)c][s7^(l-75)&] 



(5.3) 



In principle, the dependence of the Wilson coefficients Ci 
on the scale fi should cancel against the /u-dependence of 
the matrix elements of the Oi. As it is well known, the 
calculation of matrix elements such as {DsjD\Hef f\B) 
is a difficult task. The simplest evaluation can be ob- 
tained by naive factorization [2ll |. in which the — )■ 
D'^ Dsj{270Q) decay amplitude is written as: 

AiB+ ^ D°Dsji2700)) - ai^KfeKl x 

< D°{v')\bj''{l - j5)c\B{v) > X 

<i?,j(p,e)|c7^(l- 75)10 > (5.4) 

with ai = Cl ± C2/3, Gp the Fermi constant and Vij 
elements of the Cabibbo Kobayashi Maskawa mixing ma- 
trix; we use ai — 1.1 together with the values quoted in 
m for V,j. 

In the heavy quark limit, the first matrix element in 
(|5.4p can be written in terms of the Isgur-Wise function 



(5.5) 

with V and v' the four- velocities of the heavy mesons. 
The linear parametrization: 



^(v-v') = I - p^{v -v' - I) 



(5.6) 



involves the slope p^, for which various determinations 
are available: = 0.83tJj }5iE5;g^ from lattice QCD [H, 



: 1.26 ± 0.16 ± 0.11 from a Belle mea- 
The Heavy Flavour Averaging Group 



ment |2j|, / 
surement [2 

quotes — 1.23 ± 0.05 from B ^ D* decays, and 
p2 = 1.17 ± 0.18 from B D transitions ^: we use 
this last value. Computing B{B+ D°Dsj{27m)) by 
naive factorization, considering the branching fraction 
6(i:',,7(2700) D°K+) in Table HI for the two possi- 
ble assignments to Dsj(2700), and comparing the result 
to (11.211. we obtain: 



/d- = (243 ±41) MeV 
/d;^ = (181 ±30) MeV 



(5.7) 



depending on the assignements. In both cases, the lep- 
tonic constant turns out to be sizeable, similar to the 
measured decay constant: fo^ = 274±13±7 MeV 



VI. CONCLUSIONS 

We have discussed possible ways to distinguish be- 
tween two assignments for the state Dsj (2700). Our 
result is that the decay mode to D*K has very differ- 
ent branching ratios in the two possible assignments, so 
that a measurement of such a branching fraction would 
be useful to shed light on the identification of Dgj {2700). 

We have also obtained several predictions, namely the 
effective couplings governing the strong decays of D^j to 
D{D*) and a light pseudoscalar meson, which are dif- 
ferent according to the adopted interpretation. We have 
also obtained predictions for non strange partners of D^j, 
for which we also found that the decay to D*tt has the 
largest branching ratio, and predictions for the spin part- 
ners. Further reasearch on Z^sj (2700) according to these 
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TABLE IV: Branching ratios of the spin partner of Dsj for the two quantum number assignments. 





B{Di{D:^)~.D'°K+) 


B{D',{Dt2)^ D*+Ks) 




D's (J^=0-) 


(50.0 ± 0.5)% 


(23.7 ±0.2)% 


(2.6 ± 0.9)% 


Dt2 {J'' = 2-) 


(49.8 ± 0.6)% 


(23.6 ±0.2)% 


(3.1 ±1.0)% 



suggestions would be useful to complete our understand- MRTN-CT-2006-035482, "FLAVIAnet" . 
ing of the open charm meson spectrum. 
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